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ABSTRACT 



o 



, Context. After years of modest optical activity, the quasar-type blazar 4C 38.41 (B3 1633+382) experienced a large outburst in 201 1, which was 
. , detected throughout the entire electromagnetic spectrum, renewing interest in this source. 
»Cj Aims. We present the results of low-energy multifrequency monitoring by the GASP project of the WEBT consortium and collaborators, as well 
[ as those of spectropolarimetric/spectrophotometric monitoring at the Steward Observatory. We also analyse high-energy observations of the Swift 
and Fermi satellites. This combined study aims to provide insights into the source broad-band emission and variability properties. 
Methods. We assemble optical, near-infrared, millimetre, and radio light curves and investigate their features and correlations. In the optical, we 
also analyse the spectroscopic and polarimetric properties of the source. We then compare the low-energy emission behaviour with that at high 
energies. 

Results. In the optical-UV band, several results indicate that there is a contribution from a quasi-stellar-object (QSO) like emission component, 
in addition to both variable and polarised jet emission. In the optical, the source is redder- when-brighter, at least for R ^ 16. The optical spectra 

■ display broad emission lines, whose flux is constant in time. The observed degree of polarisation increases with flux and is higher in the red 

■ than the blue. The spectral energy distribution reveals a bump peaking around the U band. The unpolarised emission component is likely thermal 
' radiation from the accretion disc that dilutes the jet polarisation. We estimate its brightness to be Rqso ~ 17.85-18 and derive the intrinsic jet 

r~^. , polarisation degree. We find no clear correlation between the optical and radio light curves, while the correlation between the optical and y-ray 
flux apparently fades in time, likely because of an increasing optical to y-ray flux ratio. 
• Conclusions. As suggested for other blazars, the long-term variability of 4C 38.41 can be interpreted in terms of an inhomogeneous bent jet, 
• 1 where different emitting regions can change their alignment with respect to the line of sight, leading to variations in the Doppler factor S. Under 
the hypothesis that in the period 2008-2011 all the y-ray and optical variability on a one- week timescale were due to changes in 6, this would 
J . range between ~ 7 and ~ 21. If the variability were caused by changes in the viewing angle only, then 9 would go from ~ 2.6° to ~ 5°. Variations 
in the viewing angle would also account for the dependence of the polarisation degree on the source brightness in the framework of a shock-in-jet 
' model. 

^ ' Key words, galaxies: active - galaxies: quasars: general - galaxies: quasars: individual: 4C 38.41 - galaxies: jets 

'X: 

H . 

5r ■ 1- Introduction be explained by assuming that their emission is relativistically 

beamed, which occurs if the emitting plasma jet, produced by 

The Compton Gamma Ray Observatory (CGRO) launched in a erm assive black hole fed b an accret ion disc, is directed 

1991 revealed 27 y-ray sources on^fotirth of whlc h were iden- md ufj The larised , low . en ergy emission (from the radio 

lifted as blazars ( |Hartman et alj. 11223), i*. active galactic nu- tQ ±e optical _X- r ay band) is likely synchrotron radiation pro- 

clei (AGNs) showing the most extteme properties. They are in- duced fe relativistic electrons in the jet, while the high-energy 

deed characterised by violent activity at almost all frequencies emifjsion (from ±e x tQ ±e frequencies) is usually 

on different timescales, from long-term flux changes to inttaday inte ted as the result of i nve rse-Compton scattering of low- 

variability (IDV), by high radio to optical polarisation, and by hotonfj off ±e same relativistic electrons. Whether these 

the superluminal motions of radio knots. The observations can low . ener gy photons come from the jet itself (synchrotron self- 

* „, 7 ~ 7~TTZ ii ,ii ( , „ . cr) „ 7CD t ii u Compton, or SSC, models) or from the AGN environment (ex- 
The radio-to-optical data collected by the GASP-WEBT collabo- f _ _„ ... , . . , , , 

, ■ ^AonwcDT u- e j ■ ternal Compton, or EC, models), and in the latter case whether 

ration are stored in the GASP-WEBT archive; lor questions regarding Z 

their availability, please contact the WEBT President Massimo Villata the y come from the accretion disc, the broad line region, or the 

(villata@oato . inaf . it). dust y toms ' ls stl11 a matter of debate. 
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Blazars include flat spectrum radio quasars (FSRQs) and 
BL Lacertae objects. Frequently FSRQs exhibit quasi-stellar- 
object (QSO) like broad emission lines, as well as a blue 
and unpolarized continuum that is presumed to be the signa- 
ture of the thermal emissio n from the accretion disc, the so- 
called "big b l ue bump" (e.g. Wills et ail 



199llPian et all[l999t 



IRaiteri et all 12007b IrTAmmando et all 120091) . In contrast, BL 
Lacertae objects may have by definition at most weak lines, even 
if sometimes these so urces challenge their classification (e.g. 
IVermeulen etall 1 1995b . 

Researchers are trying to understand the structure of 
blazars and the mechanisms behind their emission and vari- 
ability by analysing the multifrequency behaviour of these ob- 
jects extended over the broadest possible energy range with 
data sampling and quality that is as h i gh as possible (see 
e.g. iMarscher et all (201 a Uorstad et all 120101: Ugudoet all 
I2.Q1 lallbh . The international collaboration known as the Whole 
Earth Blazar Telescope (WEBTfl was created in 1997 for 
this purpose, and involves tens of optical, radio, and near- 
infrared observatories. The WEBT campaigns have produced 
low-frequency light curves of extraordinary sampling, and these 
results have often been analysed in co njunction with high- 
frequency d a ta from s atellites (see e . g. IVillata et all [2 006: 
Raiteri et ail l2008bllcl lLarionov et all l2008t iBottcher et all 
2009: IVillata et all l2009bl: IRaiteri etall 120091, and references 



therein). 

Renewed interest in blazars occurred with the launch of 
the new-generation y-ray satellites Astrorivelatore Gamma a 
Immagini Leggero (AGILE) in 2007, and particularly with that 
of Fermi (formerly GLAST) in 2008, which operates in sur- 
vey mode. To acquire low-energy data to compare with the 
high-energy observations of AGILE and Fermi, in 2007 the 
WEBT started the GLAST-AGILE Support Program (GASP; 
see e.g. IVillata et all I2008L |2009ah iD'Am mando etall l2009t 
IRaiteri et all 1201 (X boilah b'Ammando etall 1201 lb . Since 
2008, a ground-based monit oring programme h as been running 
at the Steward Observatory (Smith et al., 2009), providing sup- 
port for the Fermi gamma-ray telescope. It uses the 2.3 m Bok 
and 1 .54 m Kuiper telesc opes with the SPOL spectropolarime- 
ter (ISchmidt et all Q_992 ) and provides publicly available spec- 
tropolarimetry, spectrophotometry, and calibrated broad-band 
flux measurements for about 40 blazarfl 

In this paper, we study the multifrequency behaviour of one 
blazar, namely the optically violent variable (OVV) FSRQ 4C 
38.41 (1633+382) at redshift z = 1.814, which belongs to the 
target list of both the GASP and the Steward Observatory mon- 
itoring programmes. The period of obs ervations presented here 
includes a big optical outburst in 201 1 (Rait eri et all 1201 lb ) that 
was also detected at y-ray energies by Fermi dSzosteki 1201 ll) . 
and in the X-ray and UV bands by the Swift satellite, providing 
an excellent opportunity to study the source variability over most 
of the electromagnetic spectrum. 

4C 38.41 had already been observed in y-rays by the 
Energetic Gamma Ray Experiment Telescope (EGRET) instru- 
ment onboard CGRO several times, and from these data its 
emission was found to vary significantly on a day timescale 
(Matt ox etall Il993l) . The maximum flux detected by EGRET 
above 100 MeV was (107.5 + 9.6) x lO^phcm^s" 1 in mid 
September 19910 AGILE observed 4C 38.41 several times. In 
particular, a preliminary analysis of the AGILE Gamma Ray 



Imaging Detector (GRID) data between 2009 December 1 and 
2010 November 30 yielded a y-ray flux FE>iooMev = (28 ± 
5) x 10~ 8 phcirT 2 s~'. This flux is consistent with the value re- 
ported in the Second Fermi Large Area Telescope (LAT) Catalog 
dNolan et all |2012|) . A refined analysis of the AGILE/GRID 
data, along with a search for possible short-term variability, 
will be presented in Vercellone et al. (2012, in preparation). 
Multifrequency observations of 4C 38.41 during the y-ray flares 
observed by Fermi in 2009-2010, including Very Long Baselin e 
Array (VLBA) images, was presented bv iJorstad et all d201 ll) . 
They conclude that high states at y-ray energies are due to inter- 
action between a disturbance travelling down the jet and the 43 
GHz VLBI core. 



2. Optical and near-infrared photometry 

As mentioned in the introduction, the GASP was started in 
2007 to perform long-term optical-to-radio monitoring of se- 
lected y-loud blazars during the observations of the AGILE and 
Fermi y-ray satellites. In the optical, GASP data are collected 
in the R band only. However, for this paper we also added data 
taken by the GASP- WEBT observers in the B, V, and / bands. 
Calibration of 4C 38.41 was achieved by performing differen- 
tial photometry with respect to stars A and B in the source field, 
whose s tandard magnitudes in the B, V, and R bands were de- 
rived by Villat a et all (Il997h . We obtained standard magnitudes 
in the / band from the Sloan Digital Sky Su rvey's r and i pho- 
tometr y, after applying the transformations of lChonis & Gaskelll 
(120081) . In this way, we derived / = 15.20 ± 0.05 for star A and 
/ = 15.06 + 0.06 for star B. The GASP near-infrared (near-IR) 
data are collected in the /, H, and K bands, where calibration 
is performed according to the photometric sequences obtained 
with the AZT-24 telescope at Campo Imperatorfl. 

The optical and near-IR light curves of 4C 38.41 in the 
period 2008-2011 are shown in Fig. [Q The GASP- WEBT 
observations were performed by the following observatories: 
Abastumani, Calar AltcQ, Campo Imperatore, Crimean, Galaxy 
View, Goddard (GRT), Lowell (Perkins), Lulin, Mount (Mt.) 
Maidanak, Roque (KVA and Liverpool), Rozhen, Sabadell, San 
Pedro Martir, St. Petersburg, Teide (IAC80 and TCS), Tijarafe, 
and Torino. The V and /?-band light curves are complemented 
by data taken at the Steward Observatory in the framework of 
the monitoring programme in support of the Fermi observations 
(see Sect. [3]). The figure also displays optical data acquired by 
the UVOT instrument onboard the Swift satellite (see Sect. 15. II ). 

Light curves were carefully checked and cleaned by reducing 
the data scattering through binning data close in time from the 
same telescope, when possible, and by removing clear outliers 
as well as data with errors greater than 0.2 mag. 

Substantial variability characterises the source behaviour, es- 
pecially in 2009 and 201 1. The near-IR light curves were more 
coarsely sampled than the optical ones, but show the same trend. 
In particular, they confirm the sharp brightness fall following the 
peak of the 201 1 outburst. The range of magnitude spanned over 
the whole period is 2.1-3.0 mag, depending on the band, and 
in general appears to be larger at lower frequencies. This is a 
well-known feature of FSRQs, which is ascribed to a mixture of 
variable emission from a jet and almost constant emission from 
an accretion disc. 



1 http : //www . oato . inaf.it/blazars/webt/ 

2 http : //j ames . as . arizona . edu/~psmith/Fermi/ 

3 |http:// cossc . gsf c . nasa . gov/cossc/egret/J 



4 http : //www . astro . spbu . ru/ staff /vlar/NIRlist . html 

5 Calar Alto data were acquired as part of the MAPCAT project 
http : //www. iaa. es/~iagudo/research/MAPCAT 
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Fig. 1. Optical and near-infrared light curves of 4C 38.41 in 2008-201 1 built with GASP-WEBT data (blue dots), Swift-UVOY data 
(red crosses), and data from the Steward Observatory (green plus signs). The UVOT v-band points have been shifted by -0.1 mag 
to match the ground-based data. 



2.1. Colour analysis 

The contribution of thermal radiation from the accretion disc is 
expected to lead to a redder-when-brighter spectral behaviour 
in the optical band. We verified this by performing a colour 
analysis. Figure [2] shows the source /?-band light curve and the 
B - R colour behaviour since 1995. Pre-GASP data were pro- 
vided from the Abastumani, Calar Alto, Crimean, Mt. Maidanak, 
Torino^, and St. Petersburg observatories. The historical light 
curve in the top panel shows that the source was very active in 
1995-1998 and subsequently faded to R > 16 until 2007, when it 
reached R ~ 15.5. Other flares peaking at R ~ 15.5-16 occurred 
in 2009 and at the beginning of 201 1, while during the big out- 



6 Data taken at the Torin o Observatory in 1995-1996 were published 
in lVillataetail d 19971) and lRaiteri et all d!998h . but for the present pa- 
per we reprocessed all the frames with both aperture photometry and 
Gaussian fitting of the PSF; in particular, we carefully re-analysed the 
observations of 1995 J une 27-28 that led to the detection of the optical 
outburst announced by Bos io et al 



burst of mid 201 1 the source returned to its previous brightness 
levels of 1995 and 1997 (R ^ 15). 

We derived B - R colour indices by coupling the highest- 
quality R and B data (with errors smaller than 0.05 and 0. 1 mag, 
respectively) acquired by the same telescope within 20 minutes. 
We obtained 372 indices, with a mean time separation between 
B and R exposures of about 6 minutes. The lack of a B-band light 
curve with sampling equivalent to the 7?-band one makes the B - 
R versus time plot rather discontinuous. In particular, most of 
the colour indices were derived during faint states of the source, 
as illustrated by the difference between the average magnitude 
of the whole dataset (R = 16.42) and that of the data used to 
get the colour indices (R = 16.98). The minimum and maximum 
values of the colour index are 0.39 and 1.04, respectively, and 
its average value is < B - R >- 0.64, with a standard deviation 
<x = 0.14. 

The dependence of the colour index on brightness appears 
complex, with a large dispersion in faint states. However, a gen- 
eral redder-when-brighter trend is recognisable for R ^ 16; in 
brighter states, the colour index may even decrease. This is rem- 
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iniscen t of the "saturation effect" first noticed by Villata ~et~ai] 
(120061) for another FSRQ: 3C 454.3. The shift from a redder- 
when-brighter to a bluer-when-brighter trend as the brightness 
increases would mark the transition from a quasar-like to a BL 
Lac-like spectral behaviour, the latter occurring when the syn- 
chrotron radiation dominates over the thermal emission. The dis- 
persion in the data points during faint states is apparently caused 
by at a given magnitude the source appearing to be redder during 
active periods than in more quiescent phases, as the comparison 
between the upper and middle panel of Fig.[2]suggests. However, 
we cannot rule out that at least part of this effect is due to an un- 
verifiable photometric offset between different datasets that do 
not overlap in time. 

Furthermore, assuming as a first approximation that the 
source spectrum in the optical band follows a power-law (F v oc 
v~ a ), we can derive the energy index a = [(B - R) - 0.32]/0.41, 
where we adopted Galactic extinction values of 0.048 and 0.030 
mag in the B and R bandfl respecti vely, and effective w ave- 
lengths as well as zero-mag fluxes by Bessell et aD (Il998l) . The 
range of colour index values cited above then translates into a 
range of energy index values from 0.17 to 1.76. Thus, in the 
usual spectral energy distribution (SED) representation log(vF v ) 
versus log v, the optical spectrum of 4C 38.41 is flat when B - R 
is low, i.e. in general during faint states, while it becomes steep 
for B - R > 0.73, which usually corresponds to bright levels. 

2.2. Rapid optical variability 

We analysed the /?-band light curve presented in Fig.[T]to search 
for the most noticeable intraday variability (IDV) episodes, in- 
volving changes of more than 0.3 mag in less than six hours. 
We found six nights where such rapidly occurring episodes were 
observed: on JD = 2455065, Tijarafe observations following 
St. Petersburg ones revealed a brightening of 0.32 mag in 2.82 
hours (0.113 mag/hour); on JD = 2455699, a very dense mon- 
itoring at the Mt. Maidanak Observatory showed a fading of 
0.31 mag in 4.15 hours (0.075 mag/hour); on JD = 2455702, 
observations at the KVA revealed a 0.39 mag brightening with 
respect to St. Petersburg observations about four hours before 
(0.099 mag/hour); on JD = 2455715, we derived a 0.34 mag 
change in 3.41 hours from the data of Goddard and Lowell 
(0.100 mag/hour); on JD = 2455756, intense monitoring at 
the Mt. Maidanak and Abastumani observatories revealed vari- 
ations of 0.34 mag in 4.7 hours (0.072 mag/hour); finally, on 
JD = 2455759 a very rapid brightening of 0.31 mag in 0.74 
hours was inferred by comparing Crimean with Mt. Maidanak 
data, which implies an impressive rate of 0.42 mag/hour. A care- 
ful check of the corresponding frames and photometry confirmed 
the results, so we conclude that this extreme variability episode 
is most like ly genuine. A similar ly rapid variability episode was 
reported bv lRaiteri et al.l d2008d) for another quasar-type blazar, 
3C 454.3, in January 2008 (0.40 mag/hour), and an eve n more 
extreme episode was discussed by iRaiteri et al.l (l2008al) for the 
same object, which brightened by 1.1 mag in 1.5 hours and then 
dimmed by 1.2 mag in one hour in December 2007. 

Moreover, 4C 38.41 exhibited extraordinary optical activity 
in 201 1, with flux variations of several tenths of a magnitude on 
a day timescale, as shown by Fig. [3] The most noticeable one 
was a dimming of 0.77 mag in 13.81 hours observed on JD = 
2455702-2455703, while variations ^ 1 mag were detected in 



7 From the NA SA/IPAC Extragalactic Database, 
http : / /ned .ipac.caltech. edu/ 
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Fig. 2. The fl-band light curve of 4C 38.41 in 1995-2011 (top 
panel). Blue dots indicate the data used to derive the B-R colour 
indices shown both as a function of time (middle panel) and as 
a function of magnitude (bottom panel). The black dotted and 
blue dashed lines in the top panel indicate the average magnitude 
for the whole dataset and the data used for the colour indices, 
respectively. In the middle panel, the blue dashed line marks the 
average colour index. 

52 hours on JD = 2455700-2455702 and in 72 hours on JD = 
2455805-2455808. 

3. Optical polarimetry, spectropolarimetry, and 
spectrophotometry 

Polarimetric data for the present paper were provided by the 
Steward, Crimean, Lowell, Calar Alto, and St. Petersburg ob- 
servatories. The Steward Observatory also supplied spectropo- 
larimetry and spectrophotometry data. 

Figure [4] shows the behaviour of the observed degree of op- 
tical polarisation P, as well as that of the position angle of the 
polarisation vector as a function of time. To solve the ±180°n 
ambiguity in 0, we consider only data for which P/cr P > 5 
and then require that the change across contiguous epochs be 
< 90°. Visual inspection of the figure reveals that in general pe- 
riods of high flux correspond to periods of high polarisation de- 
gree, so that there is some correlation between P and brightness, 
which is investigated further below. For the polarisation position 
angle, there appears to be a smooth rotation of about 200° in 
0.6 years, from 2008.4 to 2009.0, followed by an ever increas- 
ing rotation of about 270° in 0.6 years, from 2009.2 to 2009.8. 

Thereafter, © shows only smaller fluctuations between 100° 

and 40°in 2010 and 2011. Hence, the 2009 outburst occurred 
during a phase of noticeable change in 0, in contrast to the 201 1 
one. 
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Fig. 3. The /?-band light curve of 4C 38.41 in the most active 
period of 2011, showing noticeable intraday and interday vari- 
ability episodes. The total number of data points is 2046, 1686 
of which come from the Mt. Maidanak Observatory. 
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Fig. 4. Optical flux densities in the R band (top panel), per- 
centage of polarised flux (middle panel), and polarisation posi- 
tion angle (bottom). Data are from the following observatories: 
Steward (blue circles), Crimean (red diamonds), Lowell (black 
triangles), Calar Alto (green squares), and St. Petersburg (orange 
crosses). The + 180°n ambiguity in is solved by using only 
data for which P/cr P > 5 and then requiring that the change 
across contiguous epochs be < 90°. 
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Fig. 5. The flux in the observed frame of the C IV (blue squares) 
and C III] (red diamonds) emission lines as a function of time 
(left panel) and source brightness (right panel). During the 
period of monitoring, there is little evidence of significant 
emission-line variability in 4C 38.41. Data are derived from the 
ongoing monitoring programme at Steward Observatory. 



3.1. Emission lines 

With a redshift of z = 1.81, two prominent broad emission-lines 
fall within the 4000-7550 A spectral window that is sampled 
by the Steward Observatory spectrophotometry of 4C 38.41: 
CIV/U550 and CIII]/11909. In Fig. we plot the emission- 
line fluxes measured on 116 nights from 2008 October to 2012 
January against both time and the optical brightness of 4C 38.41 . 
The same trends for the line fluxe s are observed as seen for 
the B a imer lines in PKS 12 22+216 dSmifh et all 1201 ll) and 3C 
454.3 dRaiteri et al.Ll2008ch . That is, the data are consistent with 
the line fluxes being constant, the scatter in the measurements 
seen within an observing campaign (about a week long) be- 
ing roughly the same magnitude as during the entire monitor- 
ing period. Likewise, the substantial variations in the equivalent 
widths of the lines in 4C 38.41 (ranging from -83 to -11 A 
for C IV) are consistent with only the continuum having varied 
since 2008 August . The size of the broad-line region in quasars 
dKaspi et al.L [2000) and the added delay in variability caused by 
the redshift of the object make significant line variability impos- 
sible on timescales as short as a week, so the scatter in the data is 
a reflection of the systematic uncertainties in the spectrophoto- 
metric measurements. The average C IV line flux is measured to 

For C III], the average flux is 
The line fluxes were measured 
by fitting a single Gaussian to the line profiles. The Gaussian 
fits yield an average full width at half maximum (FWHM) of 
(4800 + 300) km s" 1 for C IV and (6500+ 1200) km s" 1 for C III]. 
The large uncertainty in the measured width of C III] is due to 
its relatively small equivalent width, especially when 4C 38.41 
is bright, and its larger FWHM relative to C IV is likely caused 
by the inclusion of A1III/11858 within the fitted line profile. 



be (3 .02+0.28)x 10~ 14 erg cirT 2 s" 1 
(1.07±0.15)x 10~ 14 erg cirT 2 s" 1 



3.2. Continuum flux and colour 

The Steward spectrophotometry can be used to examine how 
the optical continuum of 4C 38.41 varies with time and flux 
as shown in Fig. [2] but avoiding the non-variable emission fea- 
tures that fall within the standard photometric-filter bandpasses. 
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Fig. 6. The continuum colour of 4C 38.41 as a function of the 
source brightness. The continuum colour is approximated by the 
flux ratio of the two 400 A-wide bins centred on 4800 A and 
6600 A in the observed reference frame. 

Figure|6]shows the colour of 4C 38.41 as a function of brightness 
in the R band. The colour is determined by taking the ratio of the 
fluxes within two 400 A-wide bins centred on 4800 A ("blue") 
and 6600 A ("red"). Although the spectra extend to 7550 A, us- 
ing a redder measurement would include uncorrected terrestrial 
oxygen and water absorption features and would corrupt the de- 
termination of the continuum flux ratio. The trend in Fig. [6] of 
the continuum becoming redder as the object brightens confirms 
the broad-band filter results shown in Fig. [2] Despite a range 
in wavelength of only 640 A between the flux bins in the rest 
frame of the object, the correlation between colour and bright- 
ness is pronounced, and suggests that the dominant source of 
continuum light is different when 4C 38.41 is bright than when 
the blazar is faint. 

3.3. Polarisation and flux 

A strong correlation between the degree of optical linear polari- 
sation, P, and the /?-band optical brightness is evident from the 
data and shown in Fig. [7] Generally, 4C 38.41 is more highly po- 
larised when bright than when faint. Given the correlation iden- 
tified between optical continuum colour and brightness in Sect. 
13.21 the trend observed between P and brightness translates into 
a strong correlation between optical colour and degree of polar- 
isation. This correlation is shown in Fig. [8] with the continuum 
flux ratio of the 4800 A to 6600 A bins used as the measure of 
optical colour. High polarisation in 4C 38.41 is associated with 
a redder continuum. This result is further evidence that the dom- 
inant source of continuum flux is different when the object is 
bright than when it is faint, as suggested by the observed rela- 
tionship between the flux and colour. When bright, 4C 38.41 is 
generally optically redder and more highly polarised. 

In addition to the degree of polarisation, we searched for 
correlations between the polarisation position angle and optical 
flu x, colour, and P . No a pparent correlations were found. 

[jarstad et al] d201 ll) found that at least from June to 
November 2009 there was good agreement between the posi- 
tion angle of optical polarisation and that of polarisation in the 
43 GHz core, and between the general behaviour of P in the two 
bands, except for the time when, according to the VLBA images, 
a superluminal knot was passing through the core. At this time, 
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Fig. 7. The observed degree of optical polarisation plotted 
against the brightness of 4C 38.41 in the R band. Data are from 
the following observatories: Steward (blue circles), Crimean (red 
diamonds), Lowell (black triangles), Calar Alto (green squares), 
and St. Petersburg (orange crosses). Only data having P/crp > 5 
are included for clarity. 
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Fig. 8. The observed degree of polarisation plotted against the 
optical continuum colour (as defined in Pig. [6j of 4C 38.41. The 
polarisation and flux ratio measurements are simultaneous, with 
the polarisation data being derived from the median value within 
a 5000-7000 A bin from the Steward Observatory spectropo- 
larimetry. 



the optical polarization peaked, while P in the core was minimal 
and the optical position angle and the position angle in the core 
differ by ~ 90°. This implies that the core was optically thick, 
most likely due to passage of the knot, which is consistent with 
the knot kinematics. 



3.4. Wavelength dependence of the optical polarisation 

The striking correlations observed in 4C 38.41 between the opti- 
cal flux, colour, and polarisation have been noted in other strong 
emission-lin e blazars (see e.g. ISitko et ail 1 1985b ISmith et all 
1986. 120111) and have been found to be in the same sense. That 
is, as objects fade, they become bluer and are generally not as 
highly polarised. During optical flares, a redder and more highly 
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Fig. 9. Left panel: relationship between the ratio of observed po- 
larisation in the blue to red bins and colour, defined as in Fig. 
[6] Right panel: relationship between the difference between the 
polarisation position angles in the blue and red bins, and colour. 



polarised source of continuum flux appears to dominate the op- 
tical continuum. The short (daily) variability timescales and the 
power-law spectrum at optical wavelengths of the variable com- 
ponent indicate that this is the beamed synchrotron continuum 
from the relativistic jet in these sources. As a flaring object fades, 
one can wonder whether the trends seen between flux, colour, 
and polarisation are simply a reflection of the evolution of the jet, 
or if fainter emission contributions, not directly associated with 
the beamed continuum, are becoming more dominant. A key to 
understanding these trends has been investigations of how the 
linear polarisation varies with wavelength. We used 144 avail- 
able spectropolarimetric observations of 4C 38.41 obtained be- 
tween 2008 October and the end of 2012 January to examine the 
optical spectra of P, 0, and the polarised flux. 

Figure [9] shows the ratio of the fluxes in the blue to red bins 
(defined in Sect. 13.21 ) plotted against the ratio of the observed po- 
larisation in the same bins. As the blazar becomes bluer (fainter), 
the polarisation in the blue tends to be less than P measured 
in the red bin. The strength of this wavelength dependence in 
P generally increases as the object becomes bluer (fainter), al- 
though the low levels of polarisation and flux encountered re- 
sult in relatively large uncertainties in ^sooA/^eeooA- Figure [9] 
also shows the difference between the polarisation position an- 
gles determined in the two continuum bins. Most observations 
are consistent with no wavelength dependence in 0. This gen- 
eral trend of © being constant in wavelength suggests that there 
is only a single source of polarised flux from the jet that dom- 
inates at any given time. If a second non-thermal source of po- 
larised flux that has a significantly different spectral index and 
polarisation position angle becomes bright enough to compete 
with the original source, then a strong wavelength dependenc e 
in © can occur (see e.g. lHolmes etailll984l:[s"itko et all [T984). 
although this does not seem to happe n often (ISitko et allll98l 
ISmith et all Il987t iMead et all 11990ft . Much of the wavelength 
dependence may of course be hidden if the sources of polarised 
flux have similar spectral indices and/or polarisation position an- 
gles within the spectral region being observed. 

To more clearly illustrate the spectral dependence of P for 
an object as faint as 4C 38.41, we took a median of the 144 
available polarisation spectra obtained throughout the monitor- 
ing period. This was done by rotating each polarisation spectrum 
so that all of the polarised flux is in the q Stokes parameter. That 
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Fig. 10. The median flux (blue), polarised flux (red), and nor- 
malised q' and u' Stokes (green) spectra for 4C 38.41. These 
spectra were compiled from 144 spectropolarimetric observa- 
tions of the blazar. The polarised flux spectrum is derived by 
multiplying the flux spectrum by q' and displays a featureless 
power-law as expected from a synchrotron source. The polarised 
flux spectrum is scaled by a factor of 10 for clarity. See the text 
for a full description of the q' and u' spectra. Telluric absorption 
features are identified with © symbols. 



is, the observed q and u spectra were transformed (rotated) to 
q' and u', where u' averages to over the entire spectrum. The 
resulting q' spectra were then scaled to q' = 0.1 (10% polarisa- 
tion) within the 6400-6800 A bin and then median Stokes spec- 
tra were obtained. An important advantage of using q' instead 
of P is that the Stokes parameters have normal error distribu- 
tions, while P is biased since i t is a positive, definite quantity 
fe.g. lWardle & Kronbergl 1 1974ft . 

The median q' and u' spectra for 4C 38.41 are displayed in 
Fig. [10] and identify the likely cause of the general decrease in 
the polarisation toward the blue end of the spectrum. In addition 
to a decrease in q' for the continuum at A 5 6000 A, there are 
marked decreases in the polarisation at C IV and C III]. This in- 
dicates that the broad emission lines are not highly polarised, if 
at all. Indeed, the q' spectrum appears to be very similar to an 
inverted spectrum of a QSO down to the blended He II+O III] 
emission feature close to the red wing of C IV and suggests that 
the observed wavelength dependence is caused by the dilution of 
the polarised flux by an unpolarised spectrum that is similar to 
that of an optically selected QSO. The flat, featureless spectrum 
of u' supports the assertion that was made from Fig. [9] that any 
wavelength dependence in the optical polarisation position angle 
is transitory and typically not strong. 

Figure [10] also displays the median flux (F£) and polarised 
flux (q' x F^) spectra of 4C 38.41. Here, the nature of the source 
of polarised flux becomes evident. The median polarised spec- 
trum is featureless and well fit by a power-law with a spectral 
index -1.25, which is generally redder than the continuum of 
the median total flux spectrum. The spectral index of the po- 
larised continuum is identical to that of the synchrotron contin- 
uum under the assumption that the intrinsic jet polarisation (Po) 
is independent of wavelength over this spectral region. Since 
the synchrotron power-law is redder than the QSO-like emis- 
sion components (broad emission lines + blue continuum), P 
decreases into the blue. The amount of dilution and the strength 
of the wavelength dependence of the polarisation are both a 
function of the relative brightness between the unpolarised, non- 
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varying (at least on timescales measured in months to years) 
emission from the QSO component and the highly variable non- 
thermal emission from the relativistic jet and its spectral index. 
This model not only successfully explains the wavelength de- 
pendence of the observed polarisation, but also explains the cor- 
relations seen in 4C 38.41 between continuum flux, colour, and 
level of polarisation because the bluer QSO component con- 
tributes more strongly to the optical spectrum when the object is 
faint. The same general model has also been shown to apply to 
lower-redshi f t blaza rs with strong emis sion lines, su c h as 3 C 345 
(ISmith et all fl986h. PKS 1546+02 7 (ISmith et all fl99l . and 
PKS 1222+216 lSmifh etall l20lTl) . The results for 4C 38.41 
extend the evidence of the two-component nature of these types 
of blazars into the ultraviolet down to C IV/il550. 

With the polarisation and flux information in hand, it is im- 
possible to distinguish the optical spectrum of the object into 
polarised power-law synchrotron and unpolarised components 
since it is impossible to break the degeneracy between the in- 
trinsic polarisation (Pq; assumed to be constant with wavelength) 
and its brightness unless the SED of unpolarised component is 
assumed. The problem, however, is reasonably constrained as 
the unpolarised component cannot be so bright that Pq is driven 
to excessively high values. The intrinsic polarisation certainly 
must be < 100%, but P < 40-45%, which is t he observed 
maxim um for blazars (see e.g. Ilmpev et"aQ Il982l iMead et all 
1990), is likely to be the upper limit. In addition, given the longer 
timescales observed for variations in AGNs that are not blazars, 
the unpolarised component cannot be brighter than the optical 
photometric minimum observed in the light curves shown in Fig. 
□ and Fig. |2](roughly R ~ 17.8). 

Two methods have been used to break the Po-brightness de- 
generacy to estimate the true polarisation of the optical syn- 
chrotron continuum throughout the monitoring period of 4C 
38.41. First, the brightness of the unpolarised QSO component 
can be estimated from the strengths of the C IV and C III] emis- 
sion lines based on the line-to-continuum flux ra tios of a stan- 
dard t emplate spectrum for radio-quiet QSOs (e.g. lFrancis et all 
119911) . For CIV, this yields an P-band magnitude for the QSO 
component of 17.3, and for CIII], R ~ 18.0. The discrepancy 
is caused by the CIV/CIII] line ratio in 4C 38.41 differing sig- 
nificantly from the chosen template QSO, which cannot be ex- 
plained by reddening. Alternatively, if each spectropolarimetric 
observation of 4C 38.41 is fit with a combination of the QSO 
template and polarised power-law, ignoring the line fluxes and 
leaving the non-thermal spectral index, QSO brightness, and Pq 
as free parameters, we find an average R magnitude of 17.85 for 
the QSO. The median spectral index of the polarised power-law 
from these 106 fits is -1 .6, as compared to -1 .25 found from the 
median spectrum of q'. 

The estimate of the QSO brightness based on the strength 
of C IV can be ruled out, as it would make the overall P-band 
light curve difficult to explain because 4C 38.41 often becomes 
fainter. The estimates based on C III] and the model fits are con- 
sistent with the observed light curve. The intrinsic polarisation of 
the synchrotron continuum is given by Pq = P b s x P bs/(Pobs - 
PqsoX where P b s and P b s are the observed polarisation and 
flux within the R filter bandpass. In Fig. QT| the intrinsic polar- 
isation is plotted against the brightness of the blazar, assuming 
R = 17.85 for the unpolarised QSO component. Similar results 
are obtained if R = 18 is assumed, but with smaller correc- 
tions being required between the observed and intrinsic polarisa- 
tions. The correction of P for the non-beamed nuclear emission 
greatly lessens the correlation between P and optical brightness 
displayed in Fig. [7] The correlation is not completely destroyed 
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Fig. 11. The intrinsic polarisation of the optical synchrotron 
emission Pq (red filled symbols) as a function of the brightness 
of 4C 38.41 . For comparison, the observed polarisation measure- 
ments (Fig. [7j are also shown (blue empty symbols). Correcting 
P for unpolarised emission from a QSO-like component with 
R = 17.85 yields P (see text). 

as there are no cases of P $ 5% being observed when 4C 38.41 
is in a major optical outburst, but Pq approaching 20-30% can 
occur when the object is near its faintest flux levels during the 
monitoring period. 

4. Long-term observations at radio and millimetre 
wavelengths 

We analysed radio data provided to the GASP project by the 
Submillimeter Array (SMA, 230, 275, and 345 GH20) and by the 
radio telescopes of Medicina (5, 8, and 22 GHz), Metsahovi (37 
GHz), Noto (38 and 43 GHz), and UMRAO (4.8, 8.0, and 14.5 
GHz). Additional data at 86 and 230 GHz were supplied for this 
paper by the IRAM 30 m telescopq^J We also collected archival 
data, which were supplied by the SMA (230 GHz), Medicina 
(5, 8, and 22 GHz), Metsahovi (37 GHz), Noto (5, 8, 22 and 43 
GHz), and UMRAO (4.8, 8.0, and 14.5 GHz) telescopes, as well 
as 22, 37, and 87 GHz data acquired with the Metsahovi and 
Crim ean radio antennas and publ i shed in ISalonen et all d!987l) 
and iTerasranta etall d 19921 11998L 12004 12005|). The long-term 
radio light curves of 4C 38.41 at the best-sampled wavelengths 
are displayed in Fig. [12] starting from JD = 2449500. The radio 
flux densities are compared with the optical P-band ones. 

The variability maximum amplitude (F max / F m [ n ) decreases 
with decreasing frequency (with the exception of the 230 GHz 
light curve, which lacks data before 2002), suggesting that the 
corresponding jet emitting regions become larger. The shape of 
the major event in 2001-2003, clearly visible at 37, 22, and 14.5 
GHz, changes with wavelength: it shows two peaks of simi- 
lar intensity at 22 GHz, while at 37 (14.5) GHz the first (sec- 

8 These data were obtained as part of the normal m onitoring pro- 
gramme initiated by the SMA (see Gurwell et al., 2007) as well as part 
of a dedicated programme by Ann Wehrle. 

9 IRAM 30 m data were acquired as part of the POLAMI 
(Polarimetric AGN Monitoring with the IRAM-30 m-Telescope) and 
MAPI (Monitoring AGN with Polarimetry at the IRAM-30 m- 
Telescope) programme s. Data reduction was per formed following the 
procedures described in lAgudo et al U200ai2010l) . 
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Fig. 12. /?-band optical flux densities (top) compared to the radio light curves at different frequencies in 1994-201 1. Blue dots refer 
to data collected in the ambit of the GASP project, complemented by 230 GHz data from the IRAM 30 m telescope (purple crosses). 
Black plus sign symbols are historical data taken from the participating observatories' archives and from the literature (see text for 
details). 



ond) peak prevails. It is hard to say whether we are seeing 
here the evolution of a single event, or rather the effect of two 
overlapping events. Cross-correlation of the 37 and 14.5 GHz 
light curves by means of the discrete correlation function (DCF; 
Edelson & KroM 119881: iHufnagel & BregmanL 1 1992b iPetersonl 
2001 ' shows a first peak at a time lag of 120 days, and another 
peak at 360 days (see Fig. [TBI, indicating a time delay of several 
months of longer- wavelength flux variations. This is a common 
feature in blazars and suggests that lower-frequency photons are 
produced downstream in the jet from the higher-frequency ones. 



The comparison between the optical and radio-mm light 
curves (Fig. fT2l reveals a lack of general correlation: the big 
outbursts that occurred at high radio frequencies in 2001-2003 
show no optical counterpart. Moreover, the big optical outbursts 
of 1995 and 1997 do not seem to have been followed by major 
radio events on a few-month timescale. The possibility remains 
that we either missed some important optical flare or that the 
radio-optical correlation occurs on much longer timescales than 
expected for this kind of objects. While we need to wait to see 
whether a big radio outburst will follow the 201 1 optical event, 
we note that the systematically higher radio flux in 2009-201 1, 
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Fig. 13. Discrete correlation function between the 37 GHz and 
14.5 GHz light curves shown in Fig. [12] 



unusually at all wavelengths, may be related to the renewed op- 
tical activity of the past years. 

In the framework of a geometrical interpretation of blazar 
variability previously proposed for other objects (see e.g . 
Villata & Raite riLll999tlVillata et al.ll2007ll2009ttlRaiteri et all 
2009, 1201 Ol 1201 lal) . one could assume an inhomogeneous 
curved jet, where the optical and radio emitting regions were 
more misaligned in the past. In particular, a closer alignment 
of the optical region with the line of sight in 1995-1998 would 
have made the optical radiation more relativistically beamed to- 
ward us, while in 2001-2005 the radio-emitting region would 
have been better aligned, with the consequent greater Doppler 
enhancement of the radio flux. The increased activity of the past 
few years at both radio and optical frequencies would then indi- 
cate that the radio-optical misalignment has decreased and that 
both emitting regions are seen at small viewing angles. An anal- 
ysis of the variations in the viewing angle of the parsec-scale jet 
will be performed in Jorstad et al. (2012, in preparation). 



5. Swift observations 

The Swift spacecraft dGehrels et all 12004]) carries three instru - 
ments: the Burst Alert Telescope (BAT; iBarthelmv et al.l l2005). 
obser ving between 15 keV and 150 keV; the X-ray Telescope 
(XRT lBurrows et al. 1120051) . observing between 3 and lOkeV; 
and t he Ultraviolet/Optical Telescope (UVOT; Romin g et alj 
120051) . acquiring data in the 170-600 nm range. Up to 2011 
September 30, Swift targeted 69 times 4C 38.41. 

5.1. UVOT observations 

The UVOT observations were performed in the optical v, b, and 
u bands, as well as in the UV filters uvwl, uvm2, and uvw2 
dPoole et al.Ll2008l) . We reduced the data with the HEAsoft pack- 
age version 6.10, with the 20 1 1 1 30 release of the Swift/UVOTA 
CALDB. Multiple exposures in the same filter at the same epoch 
were summed with uvotimsum, and aperture photometry was 
then performed with the task uvotsource. Source counts were 
extracted from a circular region with a 5 arcsec radius, while 
background counts were estimated in a neighbouring source-free 
circular region with a 10 arcsec radius. The results are shown in 
Fig.[14](and in Fig.Q] where the v-band points have been shifted 
by -0.1 mag to match the Johnson's V-band data). The ranges of 
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Fig. 14. Optical and UV light curves of 4C 38.41 built with Swift- 
UVOT data. 



magnitude spanned in the various bands are about 2.3, 2.0, 1.6, 
1 .6, 1 .4, and 1 .4 mag going from the v to the uvw2 filter. Hence, 
the source variability decreases as the frequency increases, ex- 
tending the trend already observed in the near-IR-optical band 
in Sect. |2] 

To clarify whether there are any emission contributions in the 
UVOT energy range, we compiled SEDs for all the 36 epochs 
where good observations in all the six UVOT filters were avail- 
able. We took into a ccount the calibration updates introduced by 
iBreeveld et alJ (1201 ll) . including in particular new A e ff, count- 
rate-to-flux conversion factors, and effective areas for the UV 
filters. We calculated the (sm all) Galactic extinc tion in the vari- 
ous bands by convolving the Cardelli et al. ( 1989|) laws (setting 
R v = 3.1 and A v = 0.03'if n l) with the new effective areas, and 
used the results to obtain de-reddened flux densities. The source 
SEDs are shown in Fig. IT31 

5.2. XRT observations 

Although a few SEDs have a wavy shape, likely because of the 
data imprecision, most SEDs have a common pattern, with a 
bump peaking in the u band in the faintest states turning pro- 
gressively into a steeper and steeper curved spectrum as the 
source brightness increases. Correction for the small flux con- 
tribution of the C IV emission line (see Sect. 13. Il l to the b band 
would not change the SED shape appreciably. This again sug- 
gests that there is a thermal contribution from the accretion disc, 
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Fig. 17. The Swift-XKT spectrum of 4C 38.41 on 201 1 May 15. 
Data are fitted with a power-law with fixed Galactic absorption. 
The lower panel shows the ratio of the data to the folded model. 
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Fig. 15. Spectral energy distributions of 4C 38.41 in the optical- 
UV frequency range compiled with UVOT data. As the jet syn- 
chrotron emission becomes fainter, the contribution of the ther- 
mal radiation from the accretion disc emerges, peaking around 
the u band. 



which peaks around the u band (as expected because of the high 
source redshift) and emerges when the jet synchrotron contribu- 
tion weakens. 

We processed the XRT event files acquired in pointing mode, 
using the HEASoft package version 6.11 with the calibration 
files 20110915B We considered the observations with expo- 
sure times longer than 5 minutes, including 65 observations in 
photon-counting (PC) mode. The task xrtpipeline was run 
with standard filtering and screening criteria, in particular se- 
lecting event grades 0-12. Source counts were extracted from 
a 25 pixel circular region (~ 60 arcsec) centred on the source, 
and background counts were derived from a surrounding annu- 
lar region with radii of 100 and 150 pixels. The count rate was 
always lower than 0.5 counts/s, so no correction for pile-up was 
needed. The xrtmkarf task was used to generate ancillary re- 
sponse files (ARF), which account for different extraction re- 
gions, vignetting, and PSF corrections. The X-ray light curve 
(net count rate and de-absorbed 1 keV flux density) is shown in 
Fig. [16] and is discussed in the next section. 

We performed spectral analysis with the xspec package. 
After "grouping" the data with the task grppha in order to 
have a minimum of 20 counts per energy channel to apply the 
X 1 statistics, only ~ 20% of the spectra resulted in more than 
10 channels. We therefore used the Cash statistics, which ex- 
tends the_j^ statistics to the case o f a low number of counts 
dCashl 119791: iNousek & Shue, 1989). We applied an absorbed 



power-law model, where absorption is modelled according to 
Wilms et al. (2000) and the hydro gen column is fixed to the 
LAB0 Galactic value JV H = 0.1 1 1 x 10 21 cm -2 . We tested the re- 
sults of the Cash statistics with the goodness command, which 
performs Monte Carlo spectral simulations based on the model 
and gives the percentage of them that have fit statistics smaller 
in value than that of the data. Moreover, we compared the results 
of the Cash method with those of the^- 2 statistics, when the lat- 
ter was applicable. We obtained robust spectral fits (excluding 
the cases in the tails of the goodness distribution) in about 70% 
of cases. These cases indicate a hard spectrum, with the photon 
index F ranging from 1.31 to 1.87, with a mean value of 1.62. 
The lowest values correspond to the highest fluxes, in agreement 
with the "harder-when-b righter" trend often observed in blazars 
(e.g. lD'Ammando et al.L I20TTI and references therein). 

In particular, the data acquired on 2011 May 15 (JD = 
2455697), corresponding to the X-ray peak visible in Fig. [16] 
are best-fitted by a power-law with F = 1.31 + 0.07, 1 keV 
de-absorbed flux density fikev = (1.46 ± 0.09) //Jy, and ob- 
served!!] flux between 0.3 keV and 10 keV of (2.24 + 0.14) x 
10~ n erg cirT 2 s _1 . The result of the goodness test is 54%, which 
assures us that the observed spectrum is well-reproduced by the 
model. The corresponding spectrum and the folded model is 
shown in Fig. [171 



6. Cross-correlations from radio to y-ray 
frequencies 

Figure [16] compares 4C 38.41 light curves at different frequen- 
cies in the period 2008-2011, i.e. from the beginning of the 
Fermi y-ray observations. From top to bottom, one can see: 1) 
the daily photon flux in the 100 MeV - 300 GeV energy range; 2) 
the corresponding weekly y-ray light curv 

<Q 3) the X -ray light 



12 iKalberla etall j2005l) ; |http : //www . astro . uni-bonn . de | 

13 The unabsorbed flux is only 1% higher. 

14 The daily and weekly y-ray light curves were 
downloaded from the Fermi-LXT monitored source site 
( http : //Fermi . gsf c . nasa . gov/ ssc/data/access/lat 1, but 

their preliminary nature does not affect the results of our analysis. 

11 http://heasarc.gsfc.nasa.gov/docs/heasarc/caldb/Swift/The comparison with an independent analysis available at the ISDC 
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Fig. 16. Light curves of 4C 38.41 at different frequencies in 2008-201 1. From top to bottom: 1) the daily y-ray light curve from 
Fermi-LXT in the 100 MeV - 300 GeV range (10 _7 phcm _2 s _I ); 2) the corresponding weekly y-ray light curve with its cubic 
spline interpolation; 3) the X-ray light curve from Swift-XRT, derived as explained in Sect. 15.21 blue dots are counts s~' , while 
black squares represent 1 keV de-absorbed fluxes (/iJy), which are filled in orange when resulting from robust spectral fits; 4) the 
de-absorbed Swift-UVCfT flux density light curve in the uvwl band (mJy); 5) the de-absorbed GASP /?-band flux-density light curve 
(mJy) with the cubic spline interpolation through the 7-day binned data; 6) the millimetre radio light curve at 230 GHz (Jy). 



curve from Swift-XKT (net count rate, as well as de-absorbed 
1 keV flux densities derived as explained in Sect. 15.21 ); 4) the 
de-absorbed UV flux densities in the uvwl band (mJy); 5) the 
/?-band de-absorbed flux densities (mJy); 6) the millimetre radio 
light-curve at 230 GHz (Jy). 

A zoomed image of the most interesting periods in 2009- 
2010 and 201 1 is displayed in Figs. 1X81 and [T9l respectively. 

To make the comparison among bands easier, we highlighted 
five periods, corresponding to flaring states in y-rays. Moreover, 



http://www.isdc.unige.ch/heavens/ indeed reveals a fair 
agreement. 



we also plotted cubic spline interpolations through the weekly 
binned y-ray and optical data, in the latter case to emphasize the 
smoothed trend of the flux, beyond the fast variability. We note 
the following features: 

- In y-rays and in the optical band, where the sampling is 
dense enough, all events appear as a succession of several 
short-term flares. 

- All events (more precisely: the envelope over the substruc- 
ture) show a similar duration in y-rays and in the optical band 
(apart from event b that does not have enough optical cover- 
age). 
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Fig. 18. An enlargement of Fig. [J6]in 2009-2010. The pink ver- 
tical line marks the epoch of one of the SEDs shown in Fig. [22] 
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Fig. 19. An enlargement of Fig. [16] in 201 1. The pink vertical 
lines mark the epochs of six of the SEDs shown in Fig. [22] 



- At all frequencies, the maximum flux during the 2011 flare 
is greater than the maximum flux during the 2009 and 2010 
flares. 

- Flare a is double-peaked in both y-rays and the optical. 
However, the second y peak is higher than the first one, while 
in the optical the reverse is true. A little bump is also visible 
at 230 GHz. 

- Flare b occurred during solar conjunction so that the opti- 
cal coverage is almost null. However, the few contempo- 
raneous data do not reveal an increase in the optical flux. 
Interestingly, the millimetre light curve shows a bump simul- 
taneous with the y event. 

- The shape of flare c is similar at y-ray and optical frequen- 
cies. 

- In flare d, the most active phase in y-rays and X-rays pre- 
cedes that in the optical band and the mm flux reaches its 
maximum value. 

- A shift between the y-ray and optical smoothed-trend flux 
maxima also occurs in flare e, where furthermore the last 
optical spike has no relevant counterpart at y-ray energies. 

In conclusion, the correlation between the source emission in 
y-rays and optical (and maybe millimetre) frequencies is clear, 
because flares are observed to occur in the same periods and their 
strength increases from 2009-2010 to 201 1 at both frequencies. 
This would suggest that the jet regions where photons at these 
frequencies are produced either coincide or are very close. The 
more coarsely sampled X-ray and UV light curves indicate some 
correlation with the y-ray and optical ones, so that the X-ray and 
(the non-thermal part of) the UV radiation may also come from 
the same jet zone. 



However, the y-optical correlation is complex. The shape 
and relative brightness of the flares differ indeed considerably. 
Cross-correlation of the y-ray daily flux with the /?-band flux by 
means of the DCF yields a strong bump (DCF pea k = 0.93, see 
Fig. l20l peaking at zero time lag, but the centroid of the bump, 
which provides a more robust estimate of the time lag dPetersonL 
2001), appears shifted toward positive time-lags. Hence, the cor- 
relation between the y-ray and optical fluxes is strong, and there 
is possibly a delay in the optical variations after those in the y- 
rays, but the broadness of the DCF bump would suggest that this 
delay may be varia ble. Moreover, this result appears in contrast 
to the findings of iJorstad et al.l (1201 ll) . who found a strong cor- 
relation between the optical and y flux variations in 2009-2010, 
with the y variations lagging by 5 + 3 days. To investigate in 
greater detail the matter, we thus restrict the ca l culat ion of the 
DCF to the period investigated bv lJorstad et alJ (1201 lh . and es- 
sentially confirm their result. As shown in Fig. [21] we find a 
peak at zero time-lag, indicating that there is a strong correlation 
(DCFp ea k = 0.95) with simultaneous flux variations. However, 
the centroid of the distribution, which provides a higher time 
resolution, is r cen = -2.5 days, which means that the y-ray 
variations would follow the optical ones by 2.5 days. We de- 
termine the uncertainty in this result by performing Monte Carlo 
simulations according to the "flux redistribution/random subset 
select ion" method (FR/RSS: IPeterson et aflll998tlRaiteri et all 
2003), which tests the importance of sampling and data errors. 
Among the 1000 simulations, about 74% (k. 1 <x) of them yield 
centroids between -3 and -2 days, so we can conservatively 
conclude that the delay of the y flux variations after the optical 
ones is 2.5 + 0.5 days. If we then run the DCF for the follow- 
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Fig. 20. Discrete correlation function between the daily y-ray 
flux and 7?-band flux density light curves shown in Fig. [16] (left) 
and between the daily y-ray flux and the QSO-corrected polari- 
sation Pq (right, see Sect. [3]). 
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Fig. 21. Discrete correlation function between the daily y-ray 
flux and /?-band flux density light curves shown in Fig. [16] be- 
fore (left ) and after (r i ght) J D = 2455450 for a comparison 
with the iJorstad et alj ( 201 lb results (see the text). The blue 
squares represent the case where the light curves are stopped at 
JD = 2455763 to avoid the influence of the optical flare without 
any significant y-ray counterpart at JD ~ 2455767. 



ing period, including the 201 1 outburst, the result changes com- 
pletely (see Fig.l2TT>. The peak is still at r = days, but the corre- 
lation weakens (DCF pea k = 0.50), and the centroid indicates that 
the optical variations follow the y ones by about 7 days. This 
result is strongly affected by the presence of the optical flare 
without y-ray counterpart at JD ~ 2455767. If we indeed stop 
the light curves at JD = 2455763, the resulting DCF still shows 
a weak correlation, but now the peak and the centroid coincide 
at r = days. The reason for the lack of a strong y-optical cor- 
relation in 2011 has probably to be ascribed to the increase in 
the optical to y-ray flux ratio starting from JD ~ 2455700 and 
culminating with the last "sterile" optical flare. 

We also investigate the possible correlation between the y- 
ray flux and the QSO-corrected optical polarisation Pq derived 
in Sect. [3] The corresponding DCF, which is shown in the right 
panel of Fig. [20] does not provide evidence of any significant 
correlation. 



7. Pair-production optical depth 

In the emitting plasma jet, y-ray photons should be absorbed be- 
cause of collisions with lower-energy photons, mostly X-rays, 
which lead to pair production. Relativistic beaming of the emit- 
ted radiation can explain why we actually observe a y-ray flux. 
The minimum value of the Doppler factoiPl 5 to avoid pair pro- 
duction, i.e. t o have an optical depth r yv _ 1, was der i ved b y 
iMattox et all (119931) . and corrected by iMadejski etail d 19961) . 
under the hypothesis that the X-ray and y-ray emissions are co- 
spatial and that the emission region is spherical in the comoving 
frame 



z - vmr 2 k 2 75 t 5 x 



= [5x 10~ 4 (1 +zT za (l 



•l/(2o) 



(1) 



where a is the X-ray energy spectral index F(E) - 
Fikev(E /keV)~ a whose relationship with the photon spectral 
index r is a = r — 1, T$ is the flux variability timescale in 
units of 10 5 s, and E y is the minimum energy of the observed 
y-ray photons. By setting /175 = Hq/(15 kms^'Mpc -1 ) = 0.95, 
a = 0.31, Fi k ev = 1.46 /Jy (see Sect.ES, E y =0.1 GeV, and 
by considering a variability timescale of 12-24 hours, we obtain 
^low = 4.1-3.7. This is the lower limit to the Doppler factor af- 
fecting the radiation emitted from the jet region where the y and 
X-ray photons are produced. 

8. Broad-band SED 

Broad-band SEDs of 4C 38.41 are presented in Fig.[22]for seven 
epochs characterised by different brightness levels and good 
broad-band data coverage (see also Figs. [T8l and [T9b. They are 
obtained with simultaneous data at near-IR, optical, and higher 
frequencies, where flux variations are faster, while in the radio 
band we searched for data within a few days. Near-IR, optical, 
UV, and X-ray flux densities are corrected for Galactic extinction 
as explained in the corresponding sections. 

In the upper left inset, we show an enlargement of the near- 
IR to UV frequency range. We note that in the SEDs correspond- 
ing to JD = 2455296 and JD = 2455773, where ground and 
space optical data overlap, the space data points are lower than 
the ground-based data in the V band, especially in the fainter 
SED. This agrees with our need in Fig. Q] to shift the v-band 
UVOT data points by -0.1 mag to reconcile them with the 
ground ones and does not change the conclusion of Sect. 15.11 
that the thermal-emission contribution likely peaks around the 
u band. The shape of the low-energy bump traced by the data 
shown in Fig. [22] suggests that the synchrotron peak lies in the 
far-IR-sub-mm frequency range. 

The spectral analysis of the X-ray data acquired by Swift has 
been presented in Sect. 15.21 All the five X-ray spectra shown in 
Fig. [22] are the result of statistically valid spectral fits, whose 
spectral indices are indicated in the figure. 

To obtain the y-ray spectra, we reduced the publicly avail- 
able Fenni-LAT dato corresponding to the epochs of interest. 
The data reduction is based on the unbinned likelihood analy- 
sis with python using the ScienceTools software package ver- 
sion 

v9r23plB 

We extracted the data from a circular region of 



15 6 = [r b (l-/?cos where r b is the Lorentz factor of bulk motion 
Tb = (1 - j8 2 )~' /2 , /? is the ratio of the plasma velocity to the speed of 
light, and 9 is the viewing angle. 

16 http://Fermi.gsfc.nasa.gov/ 

17 http : //f ermi . gsf c . nasa . gov/ ssc/ data/ analysis/scitools/ 
python_tutorial . html 
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Fig. 22. Spectral energy distributions of 4C 
38.41 from the radio to the y-ray frequen- 
cies. Different colours correspond to different 
epochs, whose JD - 2450000 is indicated in the 
lower inset together with the y and X-ray spec- 
tral indices. The near-IR to UV spectral band is 
enlarged in the upper inset. Filled (empty) cir- 
cles display ground (space) data. Radio to op- 
tical ground data were collected in the frame- 
work of the GASP and collaborators; space op- 
tical, UV, and X-ray data are from Swift, y-ray 
data from Fermi. 



interest (Rol) of radius 10° centred around the location of 4C 
38.41 and included events with high photon probability by keep- 
ing only the "source class" events (event class 2). The energy 
range used in data processing is from 200 MeV to 200 GeV. We 
also filtered the data by excluding photons observed at a zenith 
angle of > 100° to reduce contamination from the Earth limb 
y-rays that result from cosmic rays interacting with the upper 
atmosphere. 

The y-ray flux and spectrum were calculated us- 
ing the current public Galactic diffuse emission model 
(gal_2yearp7v6_v0.fits) and the recommended isotropic mode{3 
for the extragalactic diffuse emission (iso_p7v6source.txt), 
which includes the residual cosmic-ray background. The instru- 
ment response function (IRFs) valid for the isotropic spectrum 
(P7SOURCE_V6) were used for data sets with front and back 
events combined to carry out the spectral analyses. 

For each of the seven epochs, we calculated the livetime 
cube, exposure map, and diffuse response of the instrument, 
and applied the algorithm as prescribed by the unbinned like- 
lihood analysis with python. We used the user contributed tool 
(make2FGLxml.py) to correctly model the background by in- 
cluding all of the sources of interest within the Rol of 4C 38.41, 
and generated the corresponding XML model file. We carried 
out the first spectral fit by setting the spectral model correspond- 
ing to 4C 38.41 to "PowerLaw2" in the file and keeping the nor- 
malisation or integral factors and photon indices of all sources 
along with the normalisation factors of the Galactic diffuse and 
isotropic components as free parameters. The XML model also 
includes sources between 10 and 15 degrees from the target, 
which can contribute to the total counts observed in the Rol ow- 
ing to the energy-dependent size of the point spread function of 
the instrument. For these additional sources, normalizations and 
index w ere kept fixed to th e values of the Second Fermi LAT 
Catalog (Nola n et al.Ll2012l) . The first fit was then used to obtain 
the final spectral fit after fixing all the free parameters of this fit 
except for those of 4C 38.41. 



18 http : // f ermi . gsf c . nasa . gov/ ssc/data/ access/lat/ 
BackgroundModels . html 



We explored the significance of the y-ray signal from the 
source using the test statistic (TS) based on the likelihood ratio 
test. To ensure a high TS value, we carried out binning using 
three different time intervals: (1) 14 days centred around JD = 
2455296, 2455640, 2455704, 2455765, and 2455773; (2) 3 days 
centred around JD = 2455697, and (3) 5 days centred around 
JD = 2455755. The TS value for all seven time intervals are 
well in excess of 100 (> 10 cr). 

We used the user contributed tool (likeSED.py) to obtain the 
spectral points and butterfly plots for each of the seven y-ray 
states. We created nine customised energy bins to ensure more 
photons per bin. The tool creates an exposure map and an event 
file for each of the energy bins. We used the final fit from our 
likelihood analysis to fit individual bins and, as before, kept pa- 
rameters of all sources fixed except for the integral factor and 
photon index of 4C 38.41. The source exhibits a slight spectral 
variation during the seven time intervals. The individual photon- 
index values for all seven epochs are displayed in Fig. [22] with 
errors indicated by their respective bow-tie patterns. 

As displayed in Fig. [22] the results of the above analysis 
reveal that the y-ray spectra are always very steep and that all 
epochs are characterised by a strong Compton dominance, as 
usually found for quasar-like blazars. This also seems to be true 
when the jet emission is faint and the thermal contribution pre- 
vails in the optical band, as shown by the JD = 2455296 SED. 
Moreover, by comparing the y-ray spectra with the X-ray ones, 
we can see that higher y-ray states correspond to higher X-ray 
ones. Finally, the data suggest that the peak of the high-energy 
bump is located in the 10 21 — 10 22 Hz frequency interval, i.e. in 
the MeV range. 

9. Discussion 

As discussed in Sect. [6] in 2009 we found a strong y-optical cor- 
relation, with a delay in the y-ray fl ux variations after th e optical 
ones of 2-3 days, in agreement with lJorstad et al.l d201 ll) . In con- 
trast, in the 201 1 outburst the correlation was weak, and the flux 
changes in the optical band apparently follow those at y-rays by 
about one week. A similar time la g of the optical flux var iations 
after the y-ray ones was found by Havashida et al. (2012) when 
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analysing the light curves of another FSRQ, 3C 279. These au- 
thors explain the time delay in terms of a decreasing ratio of the 
external radiation to the magnetic energy densities along the jet, 
which shifts the location where the optical luminosity reaches a 
maximum downstream from where the y luminosity peaks. As 
we saw in Sect. [8] 4C 38.41 also has a strong Compton domi- 
nance, which is at the basis of the Havashid a et alJ (|2012) argu- 
ment, so the same explanation might hold for this source too. 

However, our analysis suggests that the apparent optical de- 
lay may be due to an increase in the optical flux relative to that 
in y rays, even leading to an optical flare without a significant 
y-ray counterpart. A similar situation was found bv lRaiteri et alJ 
( |20 11 ah for the 2008-2009 outbursts of 3C 454.3. 

We obtained a lower limit to the Doppler f actor 6\ nvl _ 3.7- 
4.1 to avoid self-absorption of the y radiation. ISavolainen et alj 
d2010l) combined apparent jet speeds derived from high- 
resolution VLBA images from the MOJAVE project with flux 
densities at millimetre wavelengths derived from monitoring 
data acquired at the Metsahovi Radio Observatory to estimate 
the jet Doppler factors, Lorentz factors V\,, and viewing an- 
gles for a sample of 62 blazars. For 4C 38.41 they derived 
Pt mv = 29.5, 6 = 21.3, , Tb = 31.1, and 9 = 2.6°. Since according 
to IJorstad et alJ (1201 ll) the y-ray emitting regi on is close to the 
millim etre one, we can investigate whether the ISavolainen et ail 
(2010) findings agree with ours. 

If the emission is isotropic in the jet rest-frame and the in- 
trinsic spectrum follows a power-law with index a, then the 
observed flux density F v (v) = 5 n+a F' y ,(v), where n — 3 for a 
discrete, essentially point-like emission region, and n — 2 for 
a smooth, continuous jet (e.g. lUrrv & Padovanil, Il995h . Taking 
into account the uncertainty in both the emitting region structure 
and the value of a, we can assume a flux density dependence on 
6 3 . If the flux density variations that we observe are only due to 
variations in 5, then 6 = 5 max (F v /F Vimax ) 1/3 . 

Un der the hypotheses that 5 max = 21.3 from lSavolainen et alJ 
(2010) and that all the variability shown by the y-ray and optical 
lightcurves in Fig. [16] on a one-week timescale is due to varia- 
tions in the Doppler factor, we find that 5 should range between 
this maximum value and £ m ; n = 8.0 to explain the weekly y- 
ray lightcurve. In the optical band, an uncertainty comes from 
the subtraction of the QSO contribution derived in Sect. [3] from 
the /?-band flux densities. If we consider the weekly-binned R- 
band lightcurve, <5 m j n = 6.7 if Rqso = 17.85, while 6 m in = 7.5 if 
^qso = 18.0. Hence, both the y-ray and optical <S m j n are consis- 
tent with, i.e. greater than, the lower limit 6\ ovl that we obtained. 

Moreover, we can also estimate the corresponding variation 

in the viewing angle 6 = arccos{[Tb 6 - l]/[ -^r 2 -15]} under 

the assumption that all the variability is of a geomet rical nature. 
Adopting Fb = 31.1 after ISavola inen et al] (l2010l) . we obtain 
m in = 2.6° corresponding to 5 max = 21.3 for both the y-ray and 
optical emitting regions, while we derive # max = 4.8° for the y 
radiation, and # max = 5.0 or 5.3° for the optical when Rqso = 
18.0 or 17.85, respectively. 

In conclusion, the hypothesis that all the flux variability 
on a weekly timescale is due to changes in the Doppler fac- 
tor does not contradict the observations. Had we considered 
shorter variability timescales, we would have obtained lower 
(5mj n values, down to <5 m i n = 4.3 in the extreme case that we at- 
tributed all the /?-band flux density variations to changes in 6 
and set Rqso = 17.85. Even this case is still compatible with 
the estimated <5i ow . However, it is likely that the most rapid flux 
changes are due to intrinsic processes of an energetic nature, 
such as shocks propagating downstream in the jet. In contrast, 



had we considered variability timescales longer than a week, 
we would have obtained weaker constrains. Hence, we sug- 
gest that Doppler factor variations of a geometric nature pro- 
vide the most likely explanation of the long-term flux variability 
of this source, as a l ready proposed for o t her blazars (see e.g. 
IVillata et all l2009bt lLarionov et all l20Tot iRaiteri et ail 1201 lal 
and references therein). 

Furthermore, in Sect. [3] we saw that there is a clear depen- 
dence of polarisation on brightness (Fig.[7]i that cannot be com- 
pletely explained in terms of a simple dilution effect from an un- 
polarised QSO-like emission component (Fig. [TTJ. The intrinsic 
polarisation variability can be interpreted by various available 
models. Here we investigate whether the changes in the view- 
ing angle discussed above to explain the flux variations can also 
provide a plausible explanation of the variations in the polarisa- 
tion degree. If we suppose that for a statistically significant part 
of time we see shock waves propagating downstream in the jet, 
because of relativistic aberration the angle between the line of 
sight and the direction of the normal to the shock-wave front can 
be expressed as 

sin 8 

\fi = arctan( ) 

r b (cose- ^i - r b 2 ) 

and the degree of polarisation 

a + 1 (1 - rf 1 ) sin 2 ijr 
~ a + 5/3 2-(l -?7- 2 )sin 2 iA 

(e.g. iHughes eta! .1 119851) . Here rj is the degree of compression 
of the shock wave and a is the optical spectral index (see Sects. 
[2]and[3]l. The results of the model are shown in Fig. [23] for an 
optical spectral index a = 1.6 (but they are weakly dependent 
on a) and different choices of rj. They can closely reproduce the 
behaviour of the intrinsic degree of polarisation of the jet emis- 
sion Po derived in Sect. [3] as a function of brightness (here, the 
intrinsic jet flux density in the R band, which is obtained after 
correcting the observed flux for the QSO-like contribution and 
Galactic extinction). Most of the data point dispersion can be ac- 
counted for by varying the compression parameter from rj — 1.15 
to 1.7. We note that choosing Rqso = 18, which was obtained as 
a possible alternative value for the unpolarised emission compo- 
nent in Sect. [3] would shift the points toward the bottom-right, 
worsening the agreement between data and model. 

10. Conclusions 

We have presented the results of a huge observing effort on the 
FSRQ 4C 38.41, carried out by the GASP-WEBT and collab- 
orators from 2007.5 to 2011.9 and by the Steward Observatory 
from 2008.8 to 2012.1. Earlier data were also collected, so that 
the optical and radio light curves cover in total about 17 years. 
Moreover, we also analysed the UV and X-ray data acquired 
in 2007-2011 by the Swift satellite and the y-ray data taken in 
2008-201 1 by Fermi. Light curves from the near-IR to the UV 
band, spanning a factor ~ 1 1 in frequency, show a quite impres- 
sive correlation, presumably because this whole range is in the 
upper part of the synchrotron bump. In the near-IR-to-UV spec- 
tral range, the presence of a QSO-like emission contribution in 
addition to the synchrotron emission from the jet, is revealed by 
several findings: 

- the maximum flux-variation amplitude decreases with in- 
creasing frequency; 
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Fig. 23. The intrinsic polarisation of the jet emission Pq derived 
in Sect. [3] versus the jet flux density in the R band. We have 
assumed an observed brightness Rqso = 17.85 to correct for the 
QSO-like unpolarised emission contribution that affects both the 
flux and polarisation degree observed. The blue lines represent 
the results of a shock-in-jet model where the direction of the 
front wave, hence the polarisation degree, vary because of the 
variation in the viewing angle. We adopted an optical spectral 
index of 1 .6, as found in Sect. [3] Values of the shock wave degree 
of compression 77 in the 1.15—1.7 range can account for most of 
the data dispersion. 

- the optical colour shows a redder- when-brighter trend; 

- the optical spectrum reveals unpolarised broad emission- 
lines with constant flux; 

- the optical polarisation is higher in the red than in the blue; 

- the SEDs display a bump peaking around the u band in faint 
states. 

This unpolarised emission component is likely thermal emis- 
sion from the accretion disc. Thermal AGN signatures have 
apparently been foun d in several b l azars based on different 
information (see e.g . iPerlman et al.l 120081 for a review, and 
iGhisellini etafl [2009). but this is maybe the first time so much 
observational evidence has been collected for a single object. 
Furthermore, we have been able to estimate the brightness of the 
unpolarised component, /?qso ~ 17.85-18, and to correct the ob- 
served degree of polarisation for its dilution effect to obtain the 
polarisation of the jet emission. This still shows a dependence 
on brightness, which is thus an intrinsic dependence. 

The analysis of the radio light curves confirms a scenario 
where radiation of increasing wavelength is emitted in progres- 
sively larger and more external jet regions. Optical and radio 
light curves do not show a general correlation, at least on month- 
year timescales. One possible explanation, which has already 
been proposed for other blazars, is that the optical and radio 
emissions come from different jet zones that have variable orien- 
tations with respect to the line of sight. Optical and radio flares 
would thus result when the corresponding emitting regions are 
more closely aligned, with a consequent enhancement of the 
Doppler factor. 

We found a general correlation between the flux variations 
in the optical band and at y-ray energies, as the optical and y- 
ray periods of increased activity coincide. However, the shape of 
the y-ray outbursts differs from that of the optical ones, and the 
flux ratios as well as the time lags between variations at the two 



frequencies change with time. The correlation was strong during 
2009, but became weaker in 201 1, likely because of an increased 
activity in the optical band. 

We analysed broad-band SEDs of 4C 38.41 built with con- 
temporaneous data in different brightness states. A careful spec- 
tral analysis of both the X-ray data from Swift and y-ray data 
from Fermi was also performed to obtain a reliable spectral 
shape in the high-energy part of the SED. All the selected epochs 
show a strong Compton dominance, even when the jet emission 
is weak and the unpolarised component clearly emerges in the 
optical band. 

We finally discussed a geometrical interpretation of the flux 
and polarisation variability, which is able to fairly account for the 
observational data. In this view, at least the long-term flux vari- 
ations can be ascribed to changes in the Doppler boosting fac- 
tor produced by changes in the viewing angle. In particular, the 
weekly y-ray and optical light curves would imply changes in 5 
and 9 in the ranges of ~ 7-21 and ~ 2.6°-5°, respectively. When 
using these values in the framework of a shock-in-jet model, 
where the direction of the wave front, and hence the polarisa- 
tion degree, changes according to the viewing angle, we could 
also account for the trend of polarisation with brightness. 
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